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Abstract 

We study possible advantages of randomized and quantum computing over 
deterministic computing for scalar initial-value problems for ordinary differen- 
tial equations of order k. For systems of equations of the first order this ques- 
tion has been settled modulo some details in [5]. A speed-up over deterministic 
computing shown in [5] is related to the increased regularity of the solution 
with respect to that of the right-hand side function. For a scalar equation of 
order k (which can be transformed into a special system of the first order), the 
regularity of the solution is increased by k orders of magnitude. This leads 
to improved complexity bounds depending on k for linear information in the 
deterministic setting, see [8]. This may suggest that in the randomized and 
quantum settings a speed-up can also be achieved depending on k. 

We show in this paper that a speed-up dependent on k is not possible 
in the randomized and quantum settings. We establish lower complexity 
bounds, showing that the randomized and quantum complexities remain 
at the some level as for systems of the first order, no matter how large k 
is. Thus, the algorithms from [5] remain (almost) optimal, even if we re- 
strict ourselves to a subclass of systems arising from scalar equations of order k. 

Key words: fc-th order initial-value problems, randomized computing, 
quantum computing, complexity. 



1 Introduction 

In the previous paper [S] we established complexity bounds for scalar initial-value 
problems for ordinary differential equations of order k in the deterministic worst-case 
setting. We showed for instance that if the right-hand side function g depends only on 
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the solution function (not on the derivatives of the solution), then the ^-complexity 
for linear information is equal to ^(l/e) 1 ^ r+fc M , where r denotes the regularity of 
g. Hence, the order k of the equation contributes significantly to the ^-complexity. In 
this paper, we study the dependence of the e-complexity on k in the randomized and 
quantum settings. So far, the complexity in the randomized and quantum settings 
was considered for systems of equations of the first order, see [5|. For right-hand 
side functions with r continuous derivatives, the e-complexity is essentially of order 
(l/e) ly ^ r+ , where cp = in the deterministic setting, (p = 1/2 in the randomized 
setting, and if = 1 in the quantum setting (for details see Section 3). A speed-up 
for systems of equations over deterministic computing is thus by 1/2 in the exponent 
in the randomized case, and by 1 in the quantum case. Since a scalar equation of 
order k can be written as a special system of k + 1 equations of the first order, the 
upper bounds from |5] remain valid in this case. 

Intuitively, for scalar equations of order k, due to the increased regularity of the 
solution, one might expect a better speed-up over deterministic computing related 
to k. As we mentioned above, an improvement dependent on k is achieved in the 
deterministic worst-case setting by passing from the standard to linear (integral) 
information. One may hope that proper randomized or quantum approximation of 
the integrals involved in the computations will lead to algorithms with improved error 
bounds dependent on k. 

We show in this paper that such an improvement is not possible, and a speed-up 
in the randomized and quantum settings is independent of k. We establish lower com- 
plexity bounds, showing that the complexity is Q (^(1 / e) 1 ^ r+1 ^ 2 ^ in the randomized 

setting, and ^(l/e) 1 ^ r+1 ' ) j in the quantum setting, no matter how large k is. Hence, 

the algorithms defined in [S] remain almost optimal even if we restrict ourselves to 
special systems arising from scalar problems of order k. 

The paper is organized as follows. We first introduce necessary definitions in 
the three settings. Then we recall for further comparison known complexity bounds 
for initial-value problems. Main results, lower bounds on the e-complexity in the 
randomized and quantum settings, are shown in Section HI 



2 Problem definition 

We consider the complexity of a problem in the following form 

u^ k \x) = g (x, u(x), u'(x), . . . , u^ q \x)) , x€.[a,b], 
uV\a) = u{, j = 0,l,...,fc-l, U 

where < q < k, g : [a, b] x R q+1 -> R, u : [a, b] ->• R (a < b). 

For r is 1 and given positive numbers Dq, . . . , D r , we consider the class of right-hand 
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side functions g defined by 

g r = {g\geCM([a,b]x 
[a,b], yel 



x £ 



>q+l 



\d j g(x, y)\ < Dj, 
j = 0,l,...r}, 



for 



(2) 



where &>g represents all partial derivatives of order j of g. 

Instead of §Bj we can write an equivalent system of differential equations of order 1 
of the form: 
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g(u(x)), x e [a,b], (3) 



\a,u 



o 



k-l 



U„ 



(4) 



with a initial conditions 

u(a) = [c 

Then the solution u{x) of ([I]) corresponds to the function u\{x). 

Before we start analyzing the complexity of (pQ), we remind basic notions connected 
with deterministic, randomized and quantum computations. We are interested in 
finding a bounded function I = l(x) that approximates the solution of flTJ. The 
construction of I is based on a certain information about the right-hand side function 
g. In the deterministic setting, we usually consider standard information, in which 
we compute values of g or its partial derivatives at some points, or linear information 
in which we know values of linear functionals of g. 

In the randomized setting the values of g or its partial derivatives can be computed 
at randomly chosen points. In the quantum setting, information about g is gathered 
by applications of a quantum query for g. The reader is referred to [3j for a detailed 
explanation what a quantum query is. 

To get an approximate solution l(x), we use an algorithm A, which is a mapping 
from the information space into the space of bounded functions. 
In the deterministic setting, the worst-case error of an algorithm A in class Q r is 
defined by 



, worst 



(A,g r 



sup sup 

g&Q r x&[a,b~ 



\u(x) — l(x)\. 



(5) 



In the randomized and quantum settings an approximation obtained is random. Let- 
ting (Q, S, P) be a probability space, an algorithm A provides us with an approximate 
solution I", where uo G Q. 

The error of the algorithm A at g is defined by 



'(A 9) 



sup 

xd fa, 6 



\u{x)-r{x)\ 



(6) 



(we assume that e u {A ) g) is a random variable for each g e Q r ). The error of A in 
the class Q r in the randomized setting is defined by 



(A,£ r ) = sup (E(e^(A,s)) 2 ) 

9&Q r 



(7) 
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and in the quantum setting by 

e quant (A,£ r ) = e quant (A,£ r ,5) = sup inf{a | P(e w (A#) > a) < 5}. (8) 

geQ r 

The number 5 G (0, 1/2) denotes here the failure probability. It is often assumed 
to be 1/4. The success probability can then be increased by taking a median of a 
number of repetitions of an algorithm A (see [3]). 

By the cost in the deterministic, randomized and quantum settings, we mean a 
number of subroutine calls for g, which is used by an algorithm A. Thus, in the 
deterministic and randomized settings, the cost is equal to number of evaluations of g 
or its partial derivatives, while in the quantum setting it is a number of quantum query 
calls. We will denote the cost in the respective setting by cost worst (v4), cost rand (74) or 
cost quant (A). 

For any e > 0, by the e-complexity of the problem we mean the minimal cost 
sufficient to solve the problem with error no larger than e, where the minimum is 
taken over all algorithms solving the problem 

comp(</ r ,e) = min {cost(A) | e(A, Q r ) ^ e} . (9) 

To denote the complexity in the deterministic, randomized or quantum settings, we 
will use a suitable superscript: "worst", "rand" or "quant". Additionally to denote 
different types of information used in the deterministic setting, we will use indices: 
"worst-st" and " worst- lin" for standard and linear information, respectively. 



3 Upper complexity bounds 

In this section we briefly recall known complexity bounds for scalar equations of 
order k, as well as those for systems of the first order. 

In [5], Kacewicz dealt with systems of equations of the first order of the form 

z\t) =/(*(*)), tE[a,b], z{a)=r), (10) 

where / : IR d — > M. d and r] G M. d . He considered the Holder class of functions 

F r 'P = {f : R d -> R d | / G CW(R d ), ^Piy)] < D i: i = 0,l,...,r, 

\d r P{y)-d r P{z)\ <H\\y-z\y\ y,zeR d , j = 1, 2, . . . , d}, (11) 

where p G (0, 1]. 

It was shown in [5] that the e-complexity is 

1 \ l/(r+p+l/2- 7 ) N 

comp iand (^,e) = 0( (-J | (12) 



and 



\ l/(r+p+l- 7 ) N 

comp quant (^,£) = O \ ( - J j (13) 
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with an arbitrarily small positive parameter 7. (The constants in the big-0 notation 
depend on 7.) These bounds are almost optimal, i.e. they essentially match lower 
bounds on the complexity. 

It is easy to see that the bounds above with p = hold for the class J rr , where 

Fr = {f:M*->R*\feCV(R d ), l&fiiy)] < D i} i = 0,l,...,r, 

yER d , j = l,2,...,d}. (14) 

For systems ffTUl) the e-complexity in the class T T is thus equal to 

l/(r+l/2- 7 ) N 

comp rand (.F» =0 [[-) I (15) 




and 




l/(r+l- 7 ) N 



comp quant (.F» = O - . (16) 



Since equation ([T]) can be transformed into a system of the first order Q, upper 
complexity bounds ffTBT) and ffTBT) are still valid for scalar equations of order k with g G 
Q r . (Although function g has unbounded components, we can consider an equivalent 
problem with a bounded right-hand side function having bounded partial derivatives 
up to the rth order.) 

The question that we deal with in this paper is whether these bounds can be 
improved. We ask if a speed-up can be achieved dependent on k due to the increased 
regularity of the solution. In some cases in the deterministic worst-case setting such 
a speed-up dependent on k can indeed be shown. We have shown in [8] for standard 
information that 

comp worst - st (£ r ,£) = ((^J ' j , (17) 

so that there is no dependence on k in this case. However, if we allow linear infor- 
mation about right-hand side function, we can achieve a better result. The use of 
integral information leads (for q — 0) to the upper bound 



comp worst - lin (£ r ,e) = O I ( - ) I . (18) 




The complexity in this case significantly depends on k. 

Intuitively, one may expect a speed-up dependent on k by replacing integrals in 
deterministic algorithms by effective randomized or quantum approximations. In the 
next section we show lower bounds on the complexity in both settings, which indicate 
that this intuition turns out not to be correct. 
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4 Lower complexity bounds in the randomized and 
quantum settings 

We show in this section the main results of this paper. We prove the following lower 
bounds on randomized and quantum complexity of equations of order k with the 
right-hand side function belonging to class Q r . Together with (jl5p and f|T6|) they 
show that the complexity of scalar equations (CQ) is independent of k. 



1 \ l/(r+l/2) 



Theorem 4.1 Let r ^ 1. For an arbitrary k 

comp rand (£ r ,£) =n l^J , (19) 

and 

f 1 \ l/(r+l) 

comp quant (^,£) = fif-J , (20) 

•where the constants in the "Q" notation depend only on the class Q r and k. 

Before proving Theorem 14.11 we give some auxiliary results. We shall consider a 
subclass Q[ of the class Q r given by 

gi = {g:[a,b]^R\geCV([a,b\) sup \gV\x) \ ^ D j: j = 0, . . . , r}. (21) 

xe[a,b] 

This class includes functions dependent on x only, that is, problem ([T]) for g e (?[ re- 
duces to iterated integration. Since G[ <Z G r , we have that comp((?[, e) ^ comp((? r , e). 

We now show properties of auxiliary functions that will be used in the proof of 
Theorem 14.11 Let 

i 

Mx) ■= J~Tiy J( x ~ ty+'sgnUr^dt (22) 

-l 

be a perfect B-spline of degree r with r knots in (—1, 1), see e.g. [TJ. Function U r is 
the Tchebycheff polynomial of second kind, and (x — t) r + l a truncated power function: 

:=fi*-* r ' ^ (23) 

10 for x < £. 

For an interval [c, d], we define by using ip r {x) a function 

(p r ([c,d),x) := a f-Tj—) Vv+i ~~ ^3^) ' x G ML ( 24 ) 

where m,- := sup IVv+iO^)! anc ^ a := mm ^j/ m j- From the properties of perfect 

xe[-i,i] •' r 

B-spline, the function </? r ([c, d], •) belongs to the class C^ r \[c,d]), and ■ ([c, d],c) = 
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(fr([c, d], d) — for j — 0, . . . , r. In the sequel we will use the functions (p r ([c, d], ■) 
on intervals [c, d] such that d — c < 2 and [c, d] C [a, 6]. Note that <p r ([c, d], •) G £?[ for 
d — c < 2 (if we define <^ r ([c, d], •) = beyond the interval [c, d]). 

In the proof of Theorem 14.11 we will need the value of the iterated integral of 
(p r ([c,d\, ■). 



r+m ' 



(-, \ r+m « 
- J , ~ x, / (1 - t) r+m sgn U r+1 (t)dt. 
2 J (r + m)\ J 
-l 

TTien 

d Vrn-l VI 

<p r {[c, d],y )dy d Vl ... dy m ^ = (d - c) r+m ^™. (25) 




Proof. Changing variables Xi = (2yi — d — c) /(d— c), i — 0, 1, . . . , m— 1 in the left-hand 
side of ff25l) . we get 



2/m-l 2/1 

L = J J ... J cp r ([c,d),y )dy dy 1 ...dy m ^ 1 

c c c 

r+m * Xm ~ 1 X1 

= a y-~2~J J J ■■■ J A+i(x )dx dxt . . . dx m _i 



-l -l -l 

1 Xm— 1 xi 1 

'd-c\ r+m f f [l 
a i — — / / . . . / — / ixq - t) , sgnty r+ iiriatax axi . . . ax m -i- 



J — ■ J (xo — t) r + sgnU r +i(t)dt dxo dxi . . . dx r , 



-l -l -l -l 



Changing the order of integration with respect to t and xq, and computing next the 
inner integral, we get 



1 Xm-\ X 2 1 

j \ r+m r r r 1 

d — c\ / / / 1 



L = a I— — j j j . . . J — j— J (xi - t) + + sgnU r+1 (t)dt dxi . . . dx m -i- 

-l -l -l -l 

We proceed similarly for 1 to get finally 

(, \ r+m -, \ 

-r) (^/ (1 - t)r+msgn ^ (t)dt 

□ 

= (d-c) r+m p™. 

For ?i6N, let h := (b — a)/n, and aj := a + ih for i — 0, 1, . . . , n. We now consider 
the following functions 

for xeKa i+1 ]^ for - = 01) n _ L (26) 
for [a u a i+ i] 



Obviously, for h < 2 functions /j belong to the class Q[ (see the explanation after the 
definition of <p r ([c,d\, ■)). 

The following lemma shows how to express the /c-fold integrals of the functions /j. 

Lemma 4.1 Let functions /j be defined as above. Then for any k > 1 and n > 1 we 

have 



b Vk-i yi 



J / • • • / Mvo)dyo dyi... djfr-! = /i r+fc £ ^ 



^(n-i-l) fe - m (27) 



m)! 



Proof. To prove the lemma we shall show that 



x Vk-i yi 




fi(yo)dy . ..dy k - 2 dy, 



k-i 



x Vk-i yi 



for x G [a, dj) 

J J . . . J Lp r ([ai,a i+1 ],yo)dyody 1 . . .dy k -i for xG[ai,a m ] 



k 



E 



Lr+m 

y ^(x-a m ) fe - m 



(28) 



' (A; — m)! 



for x G (a i+ i, 6] 



We prove this by induction with respect to k. For k = 1, we have from Fact 14.11 with 
m = 1 that 



X 

j fMdy 







for x G [a, a*) 
^ r ([ai,a i+ i],?/) dy for x6 [a^a.+i] 



for a; G (a i+1 ,b], 



so that the statement holds true. Let us assume that fl28|) holds for fc — 1. 
We consider three cases. By the definition of fi we have 



(i) for x G [a, dj 



//...//,Wd»...fed lH = o, 



(ii) for x G [a i; a i+ i] 

fi(yo)dy ...dy k _ 2 dy k -i 



x Vk-i yi 




tf r {[a h a i+1 ],y ) dy . . . dy k _ 2 dy k -i- 



a,i 
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(iii) Let x G (a i+ i,b}. We write 
x Vk-i y\ 

fi(Vo) dy - ■ .dy k - 2 dy k -i 

a a a 

Vk-i yi 



x 



(p r ([oi, Oi+i], Vo) dy ... dy k - 2 dyk-i 



<H+i 



The last term in the equation above follows from the inductive assumption. 
Integrating the last term and using Fact 14.11 with m = k to the first term, we 
get that 

x Vk-i yi ui 

r r r ur+m 

J J ...J My) dy... dy k _ 2 dy fc _ a = V+Vf+J] (fe _ m )! ^- Q w) fc ~ m 

a a a m— 

This ends the inductive proof of (I28I) . To prove Lemma |4.1[ it is sufficient to take 
x = b in 



□ 

We are ready to prove Theorem I4.ll 



Proof of Theorem Jf.,1 We first prove a lower bound on the e-complexity in the 
quantum setting. Consider the subclass Q\ of Q r . Let (ft be any quantum algorithm 
solving problem such that e quant (0,£ r ) < e (this yields that e quant (0, Q[) < e). 

We shall prove that cost quant (0) = Q ((l/£) 1/(r+1) ) . 

Note that for any function g G Q\ the problem (CQ) reduces to the computation of 
the fc-fold integral 

fc-i • x tfe - 1 h 



u[x) = 

3=0 



^2j{x - a) j + J J ... J g{t )dt dt 1 ...dt k - 1 , xe[a,b}. (29) 



a a 



We prove the lower bound by showing that the solution of a finite number of prob- 
lems (1291) with suitable right-hand side functions leads to the solution of the mean 
value problem. 

Let n G N. We divide the interval [a, b] into 2kn parts, with points a« = a + ih, 
where i = 0,1,... , 2kn and h := (b — a)/(2kn). Let fa be the functions from ff26l) 
with n := 2kn and h := h. 
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Let Aj be arbitrary numbers with |Aj| 1 for i = 0, . . . , kn — 1, and let X be 
an injective mapping X : {0, 1, . . . , kn — 1} — > {0, 1, ... , 2kn — 1} . Then obviously 

fen— 1 
i=0 

Consider fc such mappings and numbers Cj to be specified later on, j = 
0, 1, . . . , k — 1. From Lemma 14.11 we have that 

S = ^2 c j / • • • / E ^ fxj^yMvo dyi . . . dy fe _! 

J=° a a a i=0 

fe— 1 fen— 1 fc 

= E^ E A ^ r+fc E 7^v<( 2 ^ - - 1 )*~ m - 

j=0 i=0 m=l ^ 

After changing the summation order, we get that 

fen— l fc ^ fe-i / y ( '\ j_ -\ \ k ~ m 

i=0 m=l v 7 j=0 v 7 

We split the sum indexed by m into two parts. We have 

fcn-l 1 fc-1 / y / '\ I i \ fc— 1 

fen— 1 fc fc-1 / V I '\ \ \ \ k - m 

+ e a, j> - (> - ^ ■ 

i=0 m=2 v 7 j=0 v 7 

We now show that there exist mappings Xj and numbers Cj such that for any % 

E 1 * !-4r =l < 3 °) 



j=0 

and for m = 2, .... A; 



, . . . , IU, 

fc— 1 / u, / .\ . ., \ k—m 



J=0 V 7 

We now define the functions Xj for j = 0, . . . , k — I. Let 

Xj(i) = n(k-j)+i-l, i = 0,...,kn-l. (32) 
Then obviously Xj(i) G {0, . . . , 2fcn — 1}. For such and m = 1, . . . , k, we get that 

fc-1 / / .\ . t \ fc— m fc-1 / / / .\ . .\ fe— m 

n(fc — j) + z 



i=0 v 7 3=0 V 



2fcn 



fe— 1 ■ ■ \ k—m 



3=0 



2 2k 2kn 
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Let now 

fc-i ■ \ fc-i 

w (x) = Y^cA- + ±r-x\ (33) 

3=0 V / 

be a polynomial of x, x G K. We select Cj, j — 0, . . . , k — 1, to be numbers such that 

w(x) = 1. (34) 



The existence of such numbers follows from the fact that (l34j) defines a system of linear 
equations with unknown variables Cj and with a transpose Vandermonde matrix. The 
solution (numbers c 3 ) of this equation is independent of n. From (|34p we have that 

w (m \x) = 0, (35) 

for m — 1, . . . , k — 1. Hence 

fe— 1 . \ k—m 

for m = 2, . . . , k. Taking x = i/{2kn) in (1341) and (1361) we have the desired property 

E' (1 j i \ k ~ m Jl f or m = 1 
Q C] \2 + 2k'2k^) = \0 for m = 2,...,fc-l' ( } 

for alH = 0, 1, . . . , kn — 1. Consequently, (130]) and fl3~T|) hold true, and 

fc-1 6 , V ^fcn-1 



i=0 

fcn— 1 



(A;-l)r r I fcn ^ 



3=0 { { { i=0 

fen— 1 

v ' i=0 

Denoting 

fe Vk- i 

h •= I ■ ■ ^i/^(i)(?/o)dyo • • • djfc-i, (39) 

a a a 

we have from fl38l) that 

1 fcn— 1 fe— 1 

7- £ A 4 = (Jfen)-(6-a)-( fc+ -)2^ 1 (Jfc- l)!-^ Cj /, (40) 

Kn i=0 Vr j=0 



Vk-i yi u. 



fcn— 1 
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For each j, we now use the algorithm to compute an approximation Aj to Ij 
with error at most e and probability at least 3/4. This is done with cost equal to 
cost quant (0). The algorithm <fi approximates the solution u(x) of problem (CQ), in par- 
ticular it gives us an approximation of u{b). Sience the function u(x) satisfies (1291) the 
algorithm <p can be used to approximate the /c-fold integral. Repeating the algorithm 
O(logfc) times and computing the median, we improve the probability of success to 
be at least (3/4) 1 / fc . This is achieved with cost equal to (logfc ■ cost quant (0)). 
Denote C = (b - a)- { - k+r ^ r+l (k - 1)!/^. It follows from gO]) that the mean of A;'s, 

kn—l k—1 

l/(kn) ^2 Aj, is approximated by {kn) r C Y c jAj with error 

i=0 j=0 



k-1 



kn—l 



kn 



i=0 



<ei, 



(41) 



k-1 

and probability at least 3/4, where e\ := (kn) r C \ c j\ £ - The cost of doing this is 

i=o 

O (k log k ■ cost quant (0)) . We now use lower complexity bound of Nayak and Wu 
for computing the mean of kn numbers, see [7]. They showed that the cost of 
any algorithm for computing the mean must be Q(min{kn, 1/ei}) quantum queries 
for A , • • • , Xkn-i- Taking n = 6 ((l/£:) 1 ^ r+1 ^) we have that fi(min{fcn, l/^i}) = 
Q ((l/e) 1 ^ r+1 ' ) ) , which gives us the desired lower bound. (Note that a query for 

kn—l 

A , . • • , Xkn-i is also a query for g of the form g(x) = Yl \fx{i)(%))- 

i=0 

In the randomized setting we proceed in a similar way. The difference appears 
only in the part connected with computing the mean value by randomized algorithm. 
Algorithm <f> gives us an approximation Aj of each integral Ij, for j = 0, . . . , k — 1, 
such that 

\ 1/2 

The error of computing the mean value of A's by using Aj is 

2\ 1/2 



(42) 



E 



1 kn—l k—1 

-5>-(^rc£ c A 

i=0 j'=0 



E 



k-1 



(knYCj2cj(A 3 -I 3 

j=0 



1/2 



(43) 



Using the properties of the expectation and since 

2 



k-1 
j=0 



< 



k-1 

IE 



I I 



i,j=0 



\Ai Ii I + \ Aj 



(44) 
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we get 



1/2 

, fcn-1 fc-1 x - x 

8=0 J=0 



<(A;n) r C'^ ^|c iCi |^E|A-/d 2 +E|A i -J j | 2 ^ < £l , (45) 
fc-i 

where £\ = (kn) r C \cj\e. It was shown in [6] that any approximation satisfy- 

3=0 

ing (|4"5]1 must be based on Q(min{kn, (1/ei) 2 }) subroutine calls. Hence, to prove 
desired bound on the e-complexity in the randomized setting it suffices to take 
n = ((l/e) 1 '^ 1 ' 2 '). This completes the proof of Theorem 14.11 

□ 

From bounds ffl5l) and ffl6l) and Theorem 14.11 we get the following summarizing 
corollary. 

Corollary 4.1 The complexity of initial-value problem (Q]) with an arbitrary k > 1 
satisfies: 



1X l/( r +l/2- 7 ) N 



in the randomized setting 

comp rand (^,e)=0( (^) ) (46) 

and 

/ /l \ l/(r+l/2) N 

comp rand (£^) = ft (-J ). (17) 

in the quantum setting 

l/(r+l- 7 ) N 



comp quant (£ r ,e) = CM f - 



(48) 



and 

V(H-ir 



comp quant (^ r ,e) = nfQ j , (49) 

where the positive constant 7 arbitrary small, and the constants in the big-0 nota- 
tion depend on 7. 

It follows from the Corollary that the algorithms from [5] remain optimal in a subclass 
of systems arising from scalar equations of order k, for an arbitrary k > 1. 
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5 Final remarks 



We showed in this paper some limitations of randomized and quantum computing. 
The intuition that the order k of a differential equation may help in randomized and 
quantum setting (which is the case in the deterministic setting with linear informa- 
tion) turned out not to be true. Randomized computation gives us a speed-up only 
by 1/2, and quantum computation by 1 independently of k, which is the same as in 
case of the first order equations. Thus, to optimally solve the special initial value 
problem (CQ), we can transform it to system ([3]) and apply algorithms for a general 
system of order 1 presented in [5] . 
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